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ABSTRACT 


A  state-of-the-art  survey  has  been  performed  to  determine  if  any  work 
has  been  performed  which  associates  soil  dynamics  with  mechanical 
instability  in  helicopters.  None  was  found.  No  data  pertaining  to 
the  range  of  near  surface  soil  dynamic  properties  was  found.  A  method 
of  determining  the  mechanical  stability  of  a  helicopter  in  contact 
with  the  soil  has  been  developed.  A  parametric  study  has  been  per¬ 
formed  to  examine  the  UH2  helicopter  under  conditions  of  a  range  of 
soil  properties.  It  was  found  that  there  are  ranges  of  parameters 
which  will  cause  mechanical  instability.  It  is  recommended  that  a 
program  be  initiated  to  obtain  quantitative  near  surface  dynamic 
characteristics  of  representative  soils  through  field  testing. 
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FOREWORD 


The  work  presented  in  this  report  was  performed  by  Kaman  Aerospace 
Corporation  under  Contract  N00019-73-C-0460,  for  the  Naval  Air 
Systems  Command,  Department  of  the  Navy.  The  program  was  under  the 
technical  direction  of  Mr.  R.  E.  Malatino  of  the  Naval  Air  Systems 
Command.  This  work  was  performed  by  the  author  under  the  direction 
of  Messrs  R.  Jones  and  A.  Berman  of  Kaman  Aerospace  Corporation. 
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INTRODUCTION 


Mechanical  instability  or  ground  resonance  is  a  phenomenon  associated 
with  helicopters  that  can  lead  to  the  total  destruction  of  the  aircraft. 
Therefore,  the  helicopter  must  be  designed  to  be  free  from  mechanical 
instability  in  any  phase  of  its  operation  whether  it  be  airborne  or  on 
the  ground.  Ground  resonance  is  a  well  understood  problem  and  for 
operation  from  prepared  sites,  the  present  methods  of  analysis  can  be 
used  to  insure  that  the  vehicle  is  designed  to  be  stable.  Design 
problems  arise  when  a  helicopter  is  required  to  operate  from  unprepared 
sites. 

In  a  normal  mechanical  stability  analysis  it  is  assumed  that  the 
aircraft  is  resting  on  an  infinitely  rigid  base.  This  assumption 
is  valid  for  a  prepared  site  since  the  stiffness  of  the  surface  is 
much  greater  than  the  stiffness  of  the  landing  gear.  When  the  aircraft 
is  operating  from  an  unprepared  site,  the  assumption  may  not  be  valid. 
The  problem  arises  not  from  the  ground  resonance  analysis,  which  can 
take  soil  properties  into  account,  but  from  the  fact  that  very  little 
is  known  about  soil  properties  as  applied  to  ground  resonance.  A 
considerable  amoun'  of  data  is  available  pertaining  to  the  soil 
dynamics  involved  in  the  landing  of  fixed  wing  aircraft  but  virtually 
no  research  has  been  done  in  the  rotary  wing  field.  The  concerns 
relating  to  the  landing  of  fixed  wing  aircraft  are  buch  different 
than  those  relating  to  helicopters. 
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OBJECTIVES 


The  objectives  of  this  study  were  to  determine  if  any  work  has  been 
done  which  applies  directly  to  mechanical  instability  in  helicopters 
due  to  soil  dynamics  and  to  perform  a  mechanical  instability  analysis 
which  included  soil  dynamics.  The  work  was  divided  into  the 
following  phases: 

1.  State-of-the-art  survey 

2.  Soil  Dynamic  Model 

3.  Mechanical  Instability  Analysis 

The  state-of-the-art  survey  had  three  aspects.  The  first  was  to 
determine  if  any  work  pertaining  to  the  effect  of  soil  dynamics  on 
mechanical  instability  in  helicopters  has  been  done.  The  second 
required  that  existing  soil  modeling  techniques  be  examined  to 
determine  if  they  are  compatible  with  the  present  methods  of  mechanical 
instability  analysis.  The  t'lird  aspect  consisted  of  searching  the 
appropriate  literature  in  an  effort  to  find  values  for  various  soil 
properties  which  were  needed  to  perform  the  mechanical  instability  analysis. 

The  second  phase  required  that  a  soil  dynamic  model  be  developed 
which  could  be  used  as  part  of  a  mechanical  instability  analysis. 

The  two  requirements  on  the  model  were  that  it  must  represent  the 
soil  adequately  and  that  it  must  be  compatible  with  the  present  methods 
of  mechanical  instability  analysis. 

The  third  phase  of  the  study  was  the  performance  of  a  mechanical  in¬ 
stability  analysis  which  included  soil  dynamics,  to  determine  what 
range  of  soil  property  values  are  associated  with  mechanical  instability. 


2 


STATE-OF-THE-ART  SURVEY 


A  state-of-the-a"t  survey  was  performed  to  achieve  the  following: 

1.  Determine  if  any  work  concerning  mechanical  instability 
associated  with  soil  dynamics  has  been  reported. 

2.  Determine  the  state  of  the  art  of  soil  modeling  and  its 
app’icability  to  the  present  state  of  the  art  of  mechanical 
instability  analysis. 

3.  Determine  what  soil  property  data  has  been  collected  and 
how  it  was  collected. 

Four  bibliographies  were  obtained.  These  bibliographies  contained 
citations  on  well  over  1000  books,  papers,  and  reports  pertaining 
to  soil  dynamics.  Three  of  these  bibliographies  were  compiled 
specially  for  this  study,  the  fourth  (Ref.  24)  is  a  government  report. 

From  these  bibliographies,  52  papers  and  reports  and  two  books  were 
selected  as  being  worthy  of  further  investigation.  Of  the  52  papers 
and  reports,  21  (Refs.  1-21)  were  reviewed  in  some  detail. 

There  was  no  literature  found  which  made  any  association  between  soil 
dynamics  and  mechanical  instability  in  helicopters. 

Most  of  the  reports  which  were  obtained  could  be  generally  classified 
as  falling  into  one  of  three  categories:  (a)  wave  propagation  through 
soil;  (b)  tire-soil  interaction;  (c)  soil  foundation  dynamics.  The 
works  which  dealt  with  wave  propagation  were  scanned  for  soil  property 
data  which  could  be  useful  in  the  mechanical  stability  analysis.  The 
modeling  techniques  used  for  wave  propagation  were  much  different 
than  the  type  of  models  which  were  needed  for  this  study,  and  unfortunately, 
the  data  used  for  wave  propagation  models  was  not  the  type  which  was 
required  for  the  mechanical  instability  analysis.  At  first  glance,  the 
tire-soil  interaction  reports  seemed  to  contain  exactly  the  type  of 
information  that  was  sought.  However,  further  investigation  revealed 
that  these  reports  dealt  witn  traction  and  non-elastic  deformation  and 
the  soil  models  developed  die  not  model  the  properties  of  interest. 
Therefore,  these  reports  were  not  used  for  this  in-depth  study. 

The  last  type  of  report,  soil-foundation  dynamics,  was  found  to  contain 
models  most  similar  to  the  type  of  models  which  were  required  to  perform 
the  mechanical  instability  analysis.  Elastic  theory  is  used  almost 
universally  for  soil -foundation  models  for  several  reasons.  The  theory 
of  elasticity  has  a  strong  theoretical  background,  and  therefore, 
provides  a  solid  basis  from  which  to  start  (Ref.  5).  The  methods  of 
model i-.r  which  incorporate  linear  elastic  theory  lend  themselves  to 
direct  methods  of  solution  (Refs.  22-25).  Most  soils  exhibit  linear 
properties  over  small  changes  in  strain. 
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The  two  methods  of  employing  elastic  theory  most  widely  cited  in  the 
literature  are  the  elastic  half-space  theory  and  the  lumped  parameter 
method.  The  elastic  half-space  theory  treats  the  soil  as  a  semi- 
infinite  solid  with  a  rigid  plate  resting  on  the  surface.  Figure  1 
is  a  schematic  of  a  plate  resting  on  a  layered  elastic  half  space.  The 
elastic  half-space  theory  is  particularly  useful  for  investigation  of 
the  response  of  a  body  which  is  remote  from  the  source  of  vibration, 

but  is  in  contact  with  the  soil.  This  type  of  problem  can  be  solved 

using  a  finite  element  analysis.  Elastic  half-space  theory  is  also 
useful  when  analyzing  layered  soils,  the  properties  of  which  are  known 
individually  but  not  collectively. 

The  lumped  parameter  method  is  based  on  the  elastic  half-s^ace  theory 
but  rather  than  treating  the  soil  as  a  continuum,  all  the  properties 
are  lumped  into  a  single  spring,  mass,  and  damper  system  for  each 
degree  of  freedom  (Ref.  22),  Figure  2  is  a  schematic  representation 
of  a  lumped  parameter  model.  The  lumped  parameter  method  is  particularly 
useful  because  of  its  simplicity.  Numerical  solutions  are  easy  to 

obtain  using  a  linear  lumped  parameter  model.  It  must  be  noted  that 

several  simplifying  assumptions  are  made  in  the  transition  from  the 
elastic  half-space  theory  to  the  lumped  parameter  method.  It  is  assumed 
that  the  soil  is  isotropic  and  homogeneous.  A  model  using  lumped 
parameters  is  not  necessarily  linear  but  it  has  been  demonstrated  that 
linear  models  yield  very  good  results  over  small  changes  in  strain,  and 
in  some  cases,  can  even  be  used  for  very  large  changes. 

Before  examining  the  values  of  different  soil  properties,  it  is 
appropriate  to  briefly  discuss  the  phenomenon  under  study.  As  previously 
mentioned,  modern  helicopter  design  precludes  mechanical  instability  when 
the  aircraft  is  on  a  rigid  surface  but  not  necessarily  when  the  aircraft 
is  operating  from  an  unprepared  site. 

The  primary  forces  in  a  helicopter  which  lead  to  mechanical  instability 
are  caused  by  coupled  in-plane  hub  and  blade  motions.  The  geometry  of 
all  helicopters  is  such  that  the  forces  at  the  landing  gear  are  in  all 
three  reference  planes.  Therefore,  values  for  soil  properties  must  be 
for  the  near  surface  in  both  the  vertical  and  horizontal  directions. 

The  near  surface  is  the  most  difficult  to  classify  since  it  is  made  up 
of  both  soil  and  vegetation  and  extreme  variations  in  properties  can  be 
found  over  short  distances.  In  all  of  the  literature  reviewed,  the 
values  of  soil  properties  were  for  prepared  samples  or  prepared  sites. 

The  samples  were  prepared  so  that  they  would  be  homogeneous.  The  sites 
were  prepared  by  removing  the  near  surface  to  obtain  a  uniform  surface. 
Unfortunately,  data  obtained  by  either  one  of  these  methods  is  not 
representative  of  the  near  surface  except  in  those  rare  cases  where  the 
soil  is  homogeneous  over  the  first  several  feet  of  depth  and  is  not 
covered  by  vegetation. 
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The  U.  S.  Army  Corps  of  Engineers  Waterway  Experiment  Station  in 
Vicksburg,  Mississippi  was  visited  to  determine  the  type  of  soil 
testing  currently  being  done  and  the  types  of  dynamic  soil  modeling 
most  germane  to  the  situation.  Discussions  were  held  with  the  Chief 
of  Geodynamics,  Mr.  R.  Ballard,  and  Dr.  F.  McClean.  Both  Ballard 
and  McClean  confirmed  the  results  of  our  literature  search,  i.e., 
there  is  virtually  no  data  for  soil  properties  of  the  near  surface 
and,  in  general,  the  data  which  does  exist  is  for  vertical,  rotational 
and  rocking  motions,  not  for  horizontal.  They  also  stated  that  they 
felt  that  a  linear  lumped  parameter  model  would  be  valid  for  the  near 
surface  if  data  were  available. 
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SOIL  MODEL 


In  order  to  perform  a  mechanical  stability  analysis,  a  soil  model  had 
to  be  developed  which  met  two  criteria:  valid  representation  of  the 
soil  over  small  deflections,  and  compatibility  with  existing  methods 
of  ground  resonance  analysis. 

The  ground  resonance  analysis  requires  the  impedance  (force  per  unit 
response)  of  the  supporting  structure.  The  most  direct  method  of 
modeling  the  soil  would  be  to  use  measured  impedance  data  directly, 
but  since  no  data  could  be  found,  a  linear  lumped  parameter  model  was 
used  to  calculate  the  soil  impedance. 

The  second,  and  most  difficult,  aspect  of  the  soil  model  was  determining, 
the  range  of  values  that  should  be  used  for  each  of  the  parameters  in 
the  model.  As  previously  mentioned,  there  is  a  serious  lack  of  data 
about  the  near  surface.  Even  if  some  data  were  available,  it  would  very 
likely  be  biased  because  the  bulk  of  soil  data  is  taken  from  locations 
where  a  structure  is  intended  to  be  built  and  any  site  with  a  low 
spring  rate  would  be  rejected  on  sight,  or  feel,  before  any  quantitative 
tests  were  made. 

The  damping,  spring  rate,  ana  effective  mass  of  the  soil  are  not 
measured  directly.  The  parameters  which  are  frequently  reported  are 
density,  p,  Poissons  ratio,  v,  and  shear  modulus,  G.  These  values  can 
be  used  to  compute  the  values  of  the  spring  and  damper  rate  for  a 
lumped  parameter  model  using  the  following  equations  (Ref.  22): 


0 

(1) 

3.4r  2 

ir^r  ** 

(2) 

32(1  -  v)Gr 

(7  -  8v) 

(3) 

18.4(1  -  v)  - 

"(7  -  8v7“  ro  ^ 

(4) 
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where: 


Cx  =  the  horizontal  damping  rate  of  the  soil 

Cz  =  the  vertical  damping  rate  of  the  soil 

G  =  the  shear  modulus  of  the  soil 

Kx  =  the  horizontal  spring  rate  of  the  soil 

Kz  =  the  vertical  spring  rate  of  the  soil 

r  =  the  radius  of  the  base 
o 

o  =  the  mass  density  of  the  soil 
v  =  Poissons  ratio  of  the  soil 

These  equations  were  derived  for  a  rigid  circular  base  on  an  elastic 
half-space  with  uniform  properties  throughout.  In  the  elastic  half¬ 
space  theory,  the  only  mass  used  is  the  mass  of  the  footing:  therefore, 
the  effective  mass  of  the  soil  is  considered  zero. 

As  with  most  sciences,  the  phenomenon  which  is  easiest  to  classify  and 
predict  is  studied  most  diligently.  There  is  a  wealth  of  data  per¬ 
taining  to  various  types  of  sand.  It  was  felt  that  none  of  the  data 
found  was  representative  of  the  near  surface,  but  two  soil  types  which 
were  investigated  are  used  here  as  samples. 

The  first  type,  referred  to  as  type  A  i  1  this  report;  is  made  up  of 
clay  and  silt  which  was  located  in  the  Granite  Creek  Desert,  two  i.iiles 
southwest  of  Gerlach,  Nevada  (Ref.  10).  The  average  in-place  density 
is  45  lb/ft3.  The  Poissons  ratio  and  shear  modulus  are  .47  and  1370 
lb/in2,  respectively.  Using  these  vaii:->s  *'n  Equations  (1)  -  (4)  yields 
the  following  results: 

«2  =  297,781  Ib/ft 

Cz  =  135  Ib-sec/ft 

Kx  =  206,535  lb/ft 

Cx  =  63  lb-sec/ft 

The  second  type  of  soil,  which  is  referred  to  as  type  B  in  thii  report, 
is  a  fine  sand  (Ref.  10).  Data  was  taken  at  a  site  in  the  Nebraska  Sand 
Hills  about  50  miles  southeast  of  O'Neil,  Nebraska.  The  density  of  the 
sand  is  98  lb/ft3,  the  Poissons  ratio  is  .31  and  the  measured  shear 
modulus  is  9800  lb/in?.  Again,  applying  the  data  to  Equations  (1)  - 
(4)  yield  the  following  results. 
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Kz  =  1,636,173  lb/ft 
C2  =  405  lb-sec/ft 
Kx  =  1,378,729  Ib/ft 
Cx  =  231  Ib-sec/ft 

The  data  for  both  types  A  and  B  are  for  embedded  footings  so  the 
properties  of  the  near  surface  are  not  reflected  in  the  data,  but 
this  data  was  the  best  found. 
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MECHANICAL  STABILITY  ANALYSIS 


In  order  to  establish  the  range  of  soil  parameter  values  which  are 
associated  with  mechanical  instability,  ar  analysis  was  made  which 
included  soil  properties.  Rather  than  go  directly  into  the  effects 
of  soil  parameters,  it  would  be  well  to  review  the  mechanical  in¬ 
stability  problem  as  it  is  normally  approacht  I  and  then  consider  the 
effects  of  various  soil  properties.  The  UH-2  helicopter  was  used  for 
this  study,  but  the  methods  used  apply  equally  as  well  to  other  heli¬ 
copters  with  three  or  more  blades. 

The  UH-2  helicopter  has  a  fully  articulated  four-bladed  rotor  system. 

The  fully  articulated  rotor  system  has  potentially  greater  susceptibility 
to  mechanical  instability  than  other  rotor  systems  in  use  today. 

Mechanical  instability  is  a  coupling  between  in-plane  hub  natural 

frequencies  and  blade  natural  frequencies,  and  because  in-plane  hub 
natural  frequencies  are  associated  with  the  rigid  body  modes  of  the 
helicopter  on  its  landing  gear,  the  stiffness  of  the  landing  gear  in 
both  the  vertical  and  lateral  directions  are  of  primary  importance. 

The  vertical  and  lateral  spring  rates  of  the  landing  gear  and  tire 
produce  two  modes  of  coupled  roll  and  lateral  translations  of  the 
helicopter.  The  range  of  instability  associated  with  the  low  coupled 

mode  usually  occurs  below  the  operating  range,  and  is  well  damped  due 

to  shock  absorber  strut  motion.  However,  the  range  of  mechanical  in¬ 
stability  associated  with  the  high  mode  of  coupled  motion  is  usually 
not  well  damped,  because  it  involves  lateral  deflection  of  the  gear 
and  tire,  and,  therefore,  little  damping  is  obtained  from  the  struts. 

To  insure  freedom  from  mechanical  instability  in  the  operating  range, 
either  the  lateral  stiffness  of  the  gear  and  tire  should  be  soft  enough 
so  that  the  range  of  instability  occurs  well  below  the  operating  range, 
or  the  lateral  stiffness  of  the  gear  should  be  stiff  enough  so  that  the 
range  of  mechanical  instability  is  above  the  operating  range.  In  the 
design  of  the  UH-2  the  latter,  high  stiffness,  method  was  chosen  (Ref.  25). 

Because  the  normal  approach  assumes  that  the  tires  are  on  an  infinitely 
rigid  surface,  the  addition  of  soil  properties  will  reduce  the  stiffness 
in  both  the  vertical  and  lateral  directions.  The  effect  of  reducing  the 
stiffness  will  be  that  the  natural  frequencies  will  be  lowered.  The 
obvious  problem  associated  with  reducing  the  stiffness  is  that  the 
second  natural  frequency  could  be  lowered  to  such  an  extent  that  it 
drops  within  the  operating  range.  The  first  natural  frequency  will  be 
even  further  below  the  operating  range  than  it  previously  was.  The 
amount  of  damping  required  to  permit  safe  passage  through  the  natural 
frequency  during  start-up  and  shut-down  must  be  reevaluated. 
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The  means  chosen  to  perform  a  numerical  analysis  was  to  modify  a  ground 
resonance  computer  program  to  include  the  previously  discussed  soil 
model . 

Tne  computer  program  consists  of  two  parts.  The  f1 rst  part  calculates 
the  hub  mobilities  which  are  required  by  the  stability  analysis.  The 
second  part  is  the  stability  analysis  which  is  described  in  Appendix  I 

The  first  section  of  the  program  is  capable  of  handling  a  helicopter 
with  up  to  ten  restraining  points.  Each  restraining  point  can  consist 
of  up  to  four  sets  of  springs,  masses,  and  dampers  in  series.  Each  set 
can  consist  of  up  to  three  springs,  three  masses  and  three  dampers,  one 
for  each  of  the  three  directions. 

For  the  purposes  of  this  study,  three  restraining  points  were  used,  one 
for  the  tail  wheel  and  two  for  the  main  landing  gear.  The  fuselage 
attachment  is  represented  by  a  spring,  the  shock  absorber  strut  by  a 
spring  and  damper,  the  tire  by  a  spring,  and  the  soil  by  a  spring 
damper  and  mass.  Figure  3  is  a  schematic  of  the  representation  of  one 
main  landing  gear.  The  other  main  gear  and  the  tail  gear  are  similarly 
represented.  Table  1  contains  the  values  of  springs  and  dampers  used 
to  represent  the  UH-2  helicopt  r.  This  data  is  converted  to  lateral 
and  fore  and  aft  hub  mobilities  which  are  used  in  the  stability  section. 

The  printed  output  from  the  program  consists  of  the  two  neutrally  stable 
frequencies  for  each  rotor  rpm,  the  blade  damping  required  to  make  the 
system  neutrally  stable  and  the  hub  impedance  in  each  direction.  The 
program  was  run  by  varying  the  rpm  from  2  to  310  so  that  the  full 
operating  range  plus  run-up  was  considered.  The  stability  of  the 
configuration  under  analysis  is  determined  by  comparing  the  value  of 
required  blade  damping  for  stability  to  the  actual  value  of  blade 
damping.  If  the  required  damping  exceeds  the  actual  value,  the  con¬ 
figuration  is  unstable. 

A  parametric  study  was  carried  out  to  determine  the  values  of  soil 
parameters  which  have  an  impact  on  mechanical  instability.  The  soil 
model  consists  of  springs,  dampers,  and  masses.  A  series  of  runs  were 
made  to  map  the  entire  region  of  instability  due  to  soil  spring  rate, 
damping  and  mass  in  any  combination. 
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Figure  3.  Schematic  of  Right  Man  Gear  -  Front  View 
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TABLE  1 .  UH-2  PROPERTIES 


Inertial  Properties 


Mass 

=  402.98 

Slugs 

h 

=  7167.00 

SI ug-Ft 

=  20900.00 

Slug-Ft' 

l2 

=  20000.00 

SI  ug-Ft‘ 

Dynamic  Properties 

Left  Main 

Right 

Main 

Tail  Wheel 

K 

C 

K 

C 

K 

C 

Lb/Ft 

Lb-Sec/Ft  Lb/Ft 

Lb-Sec/Ft  Lb/Ft 

Lb-Sec/Ft 

Vertical -Structure 

96000 

0 

96000 

0 

0 

Vertical-Oler 

10200 

900 

10200 

900 

350 

Vertical -Tire 

0 

144000 

0 

Rtfllil 

0 

Lateral 

24986 

0 

24986 

0 

2846 

0 

Fore  &  Aft 

0 

30000 

0 

0 

0 

Locations 

x(ft) _ y  (ft ) _ z(ft) 


CG 

0 

0 

0 

Left  Main 

-5.5 

2.98 

-6.66 

Right  Main 

5.5 

2.98 

-6.66 

Tail  Wheel 

0 

-18.94 

-6.66 

Hub 

0 

0 

6.58 

_ Blade  Properties _ _ 

Blade  Mass  -  8.95  Slugs 

Static  Moment  About  Lag  Hinge  -  76.25  Slug-ft  ~ 

Moment  of  Inertia  About  Lag  Hinge  -  1148.2  Slug-ft^ 

Distance  From  Lag  Hinge  to  CG  -  8.52  ft 

Distance  From  Lag  Hinge  to  CR  -  .6875  ft 

Static  Lag  Frequency  Due  to  Damper  Spring  Rate  -  43.073  Rad/Sec 
Number  of  Blades  -  4 
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MECHANICAL  INSTABILITY  ANALYSIS  RESULTS 


The  mechanical  instability  analysis  parametric  study  was  made  to 
determine  the  range  of  possible  soil  property  values  associated  with 
mechanical  instability.  The  range  of  soil  property  values  used  for 
the  study  was  selected  to  insure  that  the  entire  range  over  which 
instability  occurs  was  covered.  These  values  were  not  obtained  from 
measured  soil  properties,  because  realistic  data  was  not  found  in  the 
literature  search.  Therefore,  the  overlap,  if  any,  of  soil  property 
values  which  are  associated  with  mechanical  instability  and  those 
values  which  represent  soil  properties  found  in  tactical  situations 
cannot  be  determined  at  this  time. 

The  properties  of  the  Uh-2  which  were  used  for  this  study  are  presented 
in  Table  1. 

The  results  of  the  parametric  study  are  presented  in  two  different  ways. 
The  first  set  of  curves.  Figures  4  to  6,  defines  the  regions  of  in¬ 
stability.  These  figures  are  plots  of  the  soil  spring  rate  versus  soil 
damping  required  for  neutral  stability  of  the  aircraft.  The  curves  are 
shown  for  two  different  effective  soil  masses.  The  area  under  each 
curve  is  a  region  of  instability. 

The  second  type  of  plot.  Figures  7  to  18,  shows  the  amount  of  blade 
damping  required  for  neutral  stability  versus  soil  spring  rate.  These 
curves  are  shown  for  several  values  of  soil  spring  rate  and  soil 
damping.  The  value  of  X8  for  the  UH-2  is  2.58,  therefore,  any  place 
where  the  curve  exceeds  2.58  is  an  unstable  condition. 

The  curves  on  the  plot  are  labeled  K^,  «z  and  KxKyKz.  These  labels 

refer  to  the  direction  of  the  soil  spring  constant.  For  instance,  a 
curve  labeled  Kx  has  a  spring  rate  in  only  the  x  direction,  the  y  and  z 
directions  havesoil  spring  rates  which  are  infinitely  stiff. 
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Figure  6.  Soil  Damping  Required  for  Neutral  Stability 
(Soil  Spring  Rate  in  x,  y,  A  z  Direction) 
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Figure  8.  Blade  Damping  Required  for  Neutral  Stability 
Vs  Soil  Spring  Rate 
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Figure  9.  Blade  Damping  Required  for  Neutral  Stability 
Vs  Soil  Spring  Rate 
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Figure  10.  Blade  Damping  Required  for  Neutral  Stability 
Vs  Soil  Spring  Rate 
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Figure  15.  Blade  Damping  Required  for  Neutral  Stability 
Vs  Soil  Spring  Rate 
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Figure  16.  Blade  Damping  Required  for  Neutral  Stability 
Vs  Soil  Spring  Rate 
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Figure  17.  Blade  Damping  Required  for  Neutral  Stability 
Vs  Soil  Spring  Rate 


CONCLUSIONS 


Based  on  the  state-of-the-art  surve.'  and  the  mechanical  instability 

analyses,  the  following  conclusions  can  be  made. 

1.  No  work  pertaining  to  mechanical  instability  in  helicopters 
associated  with  soil  dynamics  has  been  reported. 

2.  A  linear  lumped  parameter  soil  model  is  adequate  to  represent 
the  soil  for  a  mechanical  instability  analysis. 

3.  No  data  representing  the  near  surface  in  a  manner  which  was 
applicable  to  a  dynamic  model  was  found. 

4.  There  is  a  range  of  spring-rate,  damping  rate,  and  effective 
mass  of  the  soil  which  can  cause  mechanical  instability. 

5.  No  determination  could  be  made  as  to  whether  or  not  the  spring, 
damping  and  mass  values  which  can  cause  mechanical  instability 
are  representative  of  soil  conditions  which  are  found  in  tactical 
situations. 

6.  Measured  soil  impedance  data  can  be  used  to  perform  a  mechanical 
instability  analysis  for  a  helicopter  in  contact  with  the  soil. 
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RECOMMENDATIONS 


The  mechanical  properties  of  prepared  landing  sites  are  not 
well  defined  but  the  daily  operation  of  helicopters  from  such 
sites  is  proof  enough  that  they  prov-de  stable  areas  for  heli¬ 
copter  operations.  There  is  serious  question  as  to  the 
stability  of  aircraft  landing  on  some  unprepared  sites.  The 
vertical  properties  of  soil  have  been  investigated  in  some 
depth  and  empirical  relationships  exist  to  determine  what  the 
lateral  properties  are.  Unfortunately,  most  of  the  data  reported 
is  for  soil  as  a  purely  inorganic  material.  The  near  surface  on 
which  a  helicopter  would  hand  at  an  unprepared  site  contains  both 
organic  and  inorganic  material.  There  is  virtually  no  data  for 
this  type  of  soil.  Therefore,  it  is  recommended  that  a  program 
be  initiated  to  measure  the  effects  of  organic  materials  on  near 
surface  soil  properties. 

It  is  further  recommended  that  the  program  contain  the 
following  elements: 

1.  Development  of  a  method  of  obtaining  and  analyzing 
soil  data  under  both  laboratory  and  field  conditions. 

2.  Gather  data  from  two  sites  under  two  conditions: 

a.  With  the  soil  in  its  natural  undisturbed  state 

b.  With  the  near  surface  removed. 

3.  Analyze  the  data  to  determine  the  effects  of  the 
organic  materials  in  the  near  surface. 

4.  Make  recommendations  based  on  the  conclusions  drawn 
from  3. 

It  is  estimated  that  this  effort  will  require  2800  man-hours 
and  material  costs  of  $1500.00. 
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The  "ground  resonance"  type  of  mechanical  instability 
is  generalized  to  include  a  large  class  of  supporting  structures 
for  the  rotor  hub.  Coleman's  equations  for  three  or  more  blades 
are  used.  The  fuselage,  however,  is  represented  by  the  hub 
mobility  in  each  of  two  perpendicular  directions.  This  mobility 
will  vary  with  frequency  and  may  be  either  computed  or  measured 
independently  of  the  rotor  dynamics.  The  resulting  nonlinear 
equations  are  readily  solved  at  each  frequency  of  vibration  for 
the  rotor  speed  and  blade  damping  required  for  neutral  stability. 

The  generality  and  resulting  simplicity  of  the  approach 
makes  the  method  quite  suitable  for  the  analysis  of  such  diverse 
situations  as  a  helicopter  with  brakes  off  or  a  winq  mounted  rotor- 
prop  in  flight  in  various  stages  of  transition. 


NOTATION 


a  Radial  position  of  lag  hinge 

AI  -  (Axi  +  Nx1'2 

Ar  =  1  -  X2  +  (Axr  +  A^r)/2 

2  2 

AXR'AXI'  Rea^-  ant*  imaginary  parts  of  NufPX'  NwfPY 
Nr'Ni 

b  Distance  from  lag  hinge  to  center  of  mass  of  blade 


FX'  FY 
FXo'  FYo 


V  pi 


x,  V 


V  yo 


Blade  lag  damping  rate 

Forces  acting  on  hub 

Coefficients  in  =  pXQe  1U)ft,  etc 

2  2  2 

Blade  moment  of  ihertia  about  hinge  =  mbb  (1  +  r  /b  ) 
Centering  spring  rate  of  blade 
Mass  of  blade 
Number  of  blades 

Displacement  mobility  of  hub  in  x,  y  directions 
Radius  of  gyration  of  blade  about  center  of  mass 
Coordinates  of  hub  in  fixed  system 
Complex  amplitudes  of  x,  y,  x  =  xQe  ltuft,  etc 
Complex  coordinate  of  hub,  =  x  +  iy 


2  (AXI  "  AYI* 

"  2  *axr  -  Nr* 

Complex  coordinate  of  motion  of  center  of  mass 
of  rotor  =  cv  + 


5X'  CY 


Coordinates  of  motion  of  center  of  mass  of  rotor 


CXo'  CYo 


,  etc 
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A 
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Complex  amplitudes  of  cx#  cy,  ?x  =  C^e-1^ 

-  V1 

=  XB/Wf 
=  m^ab/I 

=  Kg/! 

=  A2/cj^ 

=  nmb/2(l  +  r2/b2) 

Rotor  angular  velocity 

Oj/O)^ 

Angular  whirling  velocity  measured  in  fixed 
coordinate  system 


In  1943,  the  analysis  of  the  self-excited  mechanical 
instability  of  hinged  rotor  blades  commonly  called  ground  resonance 
was  published  by  Coleman.  This  work  was  republished  in  a  NACA 
Technical  Note  in  1957^.  The  original  analysis  was  a  significant 
technological  breakthrough  and  provided  helicopter  designers  with 
an  understanding  of,  and  a  means  of  avoiding,  this  potentially 
destructive  condition. 

2  3 

Since  1943  the  theory  was  interpreted  and  simplified  '  , 

4  5  6  7 

expanded  in  scope*  ,  charted0  and  tested  .  Some  of  the  work 
made  the  theory  easier  to  use  and  some  increased  the  applicability 
but  in  no  case  has  there  been  any  significant  change  in  the  basic 
theory. 
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The  analysis  presented  in  this  paper  also  claims  no 
improvement  in  the  theory.  The  interaction  of  the  in-plane  motion 
of  the  blades  and  the  response  of  the  hub  remains  as  the  cause  of 
the  instability.  The  original  blade  equations  given  by  Coleman 
are  used  but  the  representation  of  the  hub  dynamics  is  made  very 
general  and  at  the  same  time  very  simple.  The  resulting  equations 
are  applicable  to  almost  any  structure  whose  response  can  be  cal¬ 
culated  or  measured.  The  only  limitation  is  that  at  each  frequency 
the  response  must  be  proportional  to  the  applied  force. 

It  is  the  generality  of  application  and  ease  of  solution 
which  is  claimed  as  being  novel,  but  the  basics  of  the  original 
theory  are  unchanged. 


DERIVATION  OF  EQUATIONS 


The  equations  from  Coleman  coupling  the  motion  of  the 
rotor  center  of  gravity  and  the  motion  of  the  hub  are  presented 
for  reference  (Equation  (24)  of  Reference  1) 


(m^  +  nm^) +  Bz^  +  Ba(z^  -  itoz^)  +  Kz^  +  Amz^  +  ABzf 


+  AKzf  +  nmb^  =  0 


(1) 


_2  .  2  Bq 

nmbzf  +  2nmb((l  +  -j)  (C-^  -  2iu^1  -  to 

b  mbb 


itoC1) 


,  2a_ 

+  bCx  + 


K 


6 


mbb 


tj  =  0 


21*! 


(2) 
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where  and  are  the  deflection  of  the  hub  and  the  center  of 
gravity  of  the  blades  respectively  in  complex  fixed  coordinates. 

The  first  equation  describes  the  motion  of  the  hub.  Notice  that 
the  mass  of  the  blades  are  included.  The  last  term  of  this  equation 
is  the  only  one  containing  ^  and  can  be  thought  of  as  the  negative 
of  the  force  acting  on  the  hub  due  to  the  motion  of  the  center  of 
gravity  of  the  blades.  Because  this  hub  has  only  a  single  degree 
of  freedom  in  each  direction,  the  equation  shall  be  replaced  by  a 
more  general  representation.  The  second  equation  describes  the 
motion  of  the  center  of  gravity  of  the  rotor  of  three  or  more  blades. 
The  first  term  represents  the  coupling  with  the  hub  motion.  This 
equation  is  quite  adequate  and  shall  be  retained  as  presented. 


The  variables  in  these  equations  are  complex.  It  will 
be  convenient  to  write  the  blade  equation  in  real  coordinates  using 
zf  =  x  +  iy  and  =  £x  t  icY  where  x,  y,  cx,  represent  the 
motion  of  the  hub  and  center  of  gravity  in  the  fixed  system., 

Making  these  substitutions  and  separating  the  real  and  imaginary 
parts  of  Equation  (2)  : 


*-i i - ^  x  +  CY  +  X  c 

2(1  +  r2/b2)  X  B  X 


2(1  +  r/b2)  Y  B  y 


[u2(l  -  A1)  -  A 2]CX  +  w[2 Cy  +*BCy] 
[w2(l  -  A^  -  A2Ky  -  w  [2CX+Agcx] 


=  0 


=  0 
(3) 
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The  forces  applied  to  the  hub  are 


Fx  -  ‘“Vx 


Fy  =  -nmbcY 


(4) 


For  the  condition  of  neutral  stability,  each  of  the 
variables  may  be  written  in  the  following  forms 


-ioj,t 

x  =  x  e  f 
o 

r  =  r  e“iu)ft 
A  AO 

FX  - 


(5) 


and  similarly  for  y,  Fy  where  xQ,  CXq,  FXq,  etc.  are  complex, 

representing  relative  phases,  and  is  real,  being  the  frequency 
of  the  motion  in  fixed  coordinates.  Equations  (3)  and  (4)  now 
become : 


- % — =•  +  {1  +  +  [o2(l  -  A.)  -  A,]}c  +  o)[2i  -  =  0 

2(1  +  r/b2)  6  1  2  Xo  6  Yq 


(6) 


2(1  +  r/b) 


Vrr/  {1  +  ir6  +  [“2(l '  V  '  r21)i;Yo  -  u!2i  -  V'x  -  0 


(7) 


FX0  "  “Wxo 


F„  =  nm, 

Yo  b  f  yo 


(8) 


where  the  bars  indicate  division  by  (or  in  the  case  of  A2) 


At  this  point,  the  analysis  will  deviate  from  previous 
approaches  by  describing  the  hub  dynamics  in  a  rather  general  way. 
It  is  assumed  that  the  response  of  the  hub  to  sinusoidal  forcing 
at  any  frequency  can  be  expressed  as 


(“f)FXo 


(9) 


y0  =  V“f)F*o 

where  P^,  Py  (the  displacement  mobilities)  are  complex  and  are 
functions  of  the  forcing  frequency,  to^.  This  assumption  allows 
the  possibility  of  spring  and  damping  rates  which  are  functions 
of  frequency  as  would  be  found  in  the  usual  series  arrangement  of 
tire  spring  rate  -  gear  damper  and  spring  -  fuselage  spring  rate. 
There  is  also  no  limit  on  the  number  of  normal  modes  and  natural 
frequencies  of  the  supporting  structure  that  may  be  considered  if 
these  mobilities  are  calculated.  There  is  no  limit  on  the  anisotropy 
which  may  be  considered.  It  has  been  assumed  here  that  there  is 
no  coupling  between  the  x  and  y  directions,  i.e.,  that  a  force  in 
the  x  direction  will  produce  only  motion  in  the  x  direction.  While 
this  seems  to  be  a  usually  adequate  assumption,  it  would  be  quite 
straightforward  to  extend  the  analysis  by  including  terms  of  the 
form  pxy^wf^Fyo  t*ie  expression  for  xq  in  Equation  (9). 


These  mobilities  will  be  needed  only  at  discrete  values 
of  thus  a  tabulation  of  numerical  values  will  be  perfectly 
sufficient  for  the  method  presented.  Of  course,  the  frequencies 
must  be  close  enough  so  that  no  points  of  instability  are  missed. 

The  mobilities  may  be  either  computed  or  measured.  In  either 
case,  the  mass  of  the  blades  must  be  lumped  at  the  hub  as  previously 


mentioned  in  the  discussion  of  Equation  (1)  .  The  dynamics  of  the 
blades,  in  no  way,  however,  enter  into  the  determination  of  Px  and 
Py.  These  quantities  represent  the  response  of  the  fuselage  or 
supporting  structure  as  an  independent  body. 

Typical  computational  procedures  would  consist  of 
analyzing  a  rigid  fuselage  on  the  landing  gear,  or  a  finite  element 
analysis  of  the  fuselage,  or  a  bending  analysis  of  a  wing  structure 
supporting  a  rotor-prop  in  various  positions.  Any  of  these  are 
reasonably  standard  procedures  having  no  direct  connection  with 
the  problem  at  hand  except  as  a  source  of  certain  required  data. 


When  the  structure  in  question  actually  exists,  a  rel¬ 
atively  straightforward  test  will  supply  the  needed  data.  The 
measured  response  at  the  hub  due  to  an  excitation  at  the  hub  as  a 
functin  of  frequency  is  required.  This  test  should  be  performed 
in  each  direction  and  the  complex  amplitudes  (or  amplitude  and 
phase)  should  be  recorded. 


Now  using  Equations  (8)  and  (9),  xq  and  yQ  may  be  expressed 
as  functions  of  cXq,  as  follows: 


x0  =  mVfPX(“f)?X0 
yc  =  myo*PY(uf)Cyo 


(10) 


This  may  now  be  put  into  Equations  (6)  and  (7)  to 


0 


{  [a)  (1  -  A^)  +  (1  -  A 2"*'  Ay^  ^  +  *  +  Axi^Xo  _u)  ^6  ”  2^  **Yo 

(ID 

w(Xg  -  2i)?Xo  +  {[w2(l  -  A:)  +  (1  -  I2  +  A^H  +  i(Ig  +  =  0 

(12) 


where 


rnnb 


2(1  +  r/b) 


2  2.  “fPX 


A4“fPX  =  AXR  +  lAXI 


(13) 


and  similarly  for  Py.  Note  that  A^  is  related  to  Coleman's  A^ 
through  the  effective  mass  (A^  =  A^/mg).  For  convenience  and  to 
better  represent  the  possible  anisotropy,  the  following  parameters 
are  defined. 


ar  1  “  ^2  +  2  (axr  +  Nr5 


An 


2(AXi  +  Ayi) 


2(AXR  -  Ayj^ 

2^AXI  ' 


Equations  (11)  and  (12)  then  become: 


{  [<u2  (1  -  Ax>  +  Ar  +  Ar]  +  i(Ag  +  Ax  +  Aj)}^  -  wUg  -  2i)Cyo  -  0 

(14) 

w(XB  -  2i)cXo  +  ([u2(l  -  +  Ar  -  Ar]  +  i(Ag  +  Aj  -  Ax)  >Cyo  -  0 

(15) 

The  determinant  of  the  coefficients  of  Equations  (14)  and  (15)  must 
be  zero  for  a  non-trivial  solution  to  exist.  This  is  written 

9 
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{  [<*>2  (1  -  Ax)  +  Ar  +  i(Xg  +  Aj)]  +  (Ar  +  iAj)} 

*  {  [w2  (1  -  Ax)  +  Ar  +  i  (Ag  +  Aj)]  -  (Ar  +  iAI)}  +  w2(If  -  2i) 2  =  0 

or 

[u2(l  -  Ax)  +  Ar  +  i(I6  +  Aj)]2  =  (Ar  +  iAj) 2  +  u2(2  +  iX6)2  (16) 

In  this  equation,  all  the  parameters  are  real  and  thus  the  real  and 

imaginary  parts  of  the  equation  must  separately  equal  zero.  Equa¬ 
tion  (16)  is  the  general  equation  describing  the  steady  state, 
neutraly  stable  oscillation  at  a  frequency  of  For  each  value 

of  u)^,  the  equation  may  be  solved  for  Ag  and  u)  as  will  be  described 
below. 

SPECIAL  CASE:  ISOTROPIC  SUPPORTS 


Before  proceeding  with  the  solutions  of  the  general 
equations,  it  will  be  instructive  to  examine  the  classical  but 
unusual  conditions  of  isotropy.  When  P^  =  Py,  AR  =  Aj  *  0  and 
Equation  (16)  can  be  written 

[w2(l  -  Ax)  +  Ar  +  i(A6  +  Aj)]2  =  a2  (2  +  iAg)2 


Taking  the  square  root  of  both  sides  and  separating  the  real  and 
imaginary  parts  yields 


w2(l  -  A^)  -  2w  +  Ar  =  0 

Ag(l  -  w)  +  Aj  =  0 


(17; 
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At  each  frequency  of  oscillation,  ,  AR  and  A^  are 
known  and  the  first  equation  may  be  solved  for  <U  and  then  the 
second  may  be  solved  for  Tg .  These  are  the  rotor  speed  and  the 
blade  damping  required  for  a  neutrally  stable  oscillation  at  a 
frequency 

In  the  appendix,  it  is  shown  that  for  single  degree-of- 
freedom  isotropy,  these  equations  reduce  to  forms  given  by  Coleman 
(Ref.  1)  and  Warming  (Ref.  3). 

GENERAL  CASE:  ANISOTROPIC  HUB 


For  the  general  case.  Equation  (16)  does  not  reduce  to 
so  simple  a  form  as  for  the  case  of  isotropy.  This  equation,  when 
expanded  and  when  the  real  and  imaginary  parts  are  separated,  yields 
two  equations  as  follows: 

w4(l  -  Ax)2  +  u2(Ib2  +  2(1  -  Ax)Ar  -  4}  +  {(Xg  +  Ax)2  +  Ar2 

-  (AR2  -  Aj2)}  =  0  (18) 

w2{r6(l  +  Ax)  +  Aj(l  -  Ax)}  +  (XB  +  Aj)Ar  -  ArAj  =  0  (19) 

These  equations  may  be  readily  solved  by  the  use  of  the 
Newton-Raphson  method  (See  Reference  (8),  for  example).  It  has 
been  found  that  the  solutions  for  the  isotropic  case  given  by 
Equation  (17)  are  a  good  starting  point  for  the  iteration. 
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The  results  can  be  plotted  as  in  Coleman  as  w ^  vs  to. 
A  more  useful  form,  however,  is  a  plot  of  Ag  / s  w,  showing  the 
blade  damping  required  for  stability  at  each  rotor  speed. 

ILLUSTRATIVE  EXAMPLE 


As  an  illustration  of  the  kinds  of  results  which  may 
be  obtained  using  these  procedures,  the  example  shown  here  is 
presented.  The  data  is  representative  of  a  helicopter  on  the 
ground  with  brakes  off.  Figure  1  illustrates  the  amplitudes  of 
Px  and  Py  as  a  function  of  forcing  frequency  for  two  cases,  i.e., 
with  and  without  damping  in  the  landing  gear. 

Figure  2  shows  u  as  a  function  of  for  the  two  cases. 
For  the  condition  of  no  damping,  equations  (18)  and  (19)  either 
have  no  solution  when  there  is  an  instability  or  have  a  solution 
with  Ag  •=  0  when  the  system  is  stable.  This  situation  is  well 

3 

known  .  When  there  is  no  hub  damping,  the  entire  system  is  either 
absolutely  stable  or  unstable  and  blade  damping  alone  cannot  change 
the  range  of  instability. 

Figure  3  illustrates  the  blade  damping  required  as  a 

function  of  rotor  speed.  The  ranges  of  instability  for  any  value 

of  A0  are  easily  read  from  the  figure.  For  conditions  of  lesser 
p 

gear  damping,  there  will  be  situations  where  there  will  be  no 
solution  to  the  equations.  This  indicates  a  condition  of  absolute 
instability  and  only  additional  damping  in  the  gear  can  remove 
the  instability. 
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CONCLUSIONS 


A  procedure  has  been  presented  for  applying  the  theory 
of  mechanical  instability  to  hinged  rotors  on  very  general  supports. 
The  response  of  the  hub  may  be  computed  or  measured  independently 
of  the  rotor  dynamics. 

The  resulting  procedure  is  quite  easy  to  use  and  may  be 
applied  to  almost  any  structure  supporting  a  rotor  of  three  or  more 
hinged  blades. 
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APPENDIX 


THE  ISOTROPIC  SINGLE-DEGREE-OF-FREEDOM  HUB 


The  single-degree-of-freedom  hub  on  isotropic  supports 
has  been  analyzed  by  Coleman*  and  Warming3  and  many  others.  For 


this  condition,  it  may  be  shown  that 


1  +  +  A 


—2 

wr  -  1 


2  '  “3  [ZT?  -  l]2  +  I2 


(A-l) 


3  7=5 


[Z£  -  1]2+  T2f 


(A-2) 


where  u>r  and  are  the  natural  frequency  and  damping  ratio  of  the 
hub  and  the  bars  indicate  division  by  Substituting  Equations 
(A-l)  and  (A-2)  into  Equation  (17)  results  in 


-2 


u)2(l  -  Ax) 


-  2w  +  {1  - 


l2  +  a3 


(i) 


r  -  1 


[u >2  -  l]2  +  X2 


}  =  0 


(A-3) 


6  u  -  1  [u>2  -  1]  +  X3 


(A- 4) 


At  the  "center  of  instability"  where  w  =  1,  Equation  (A-4) 
gives 


xBxf  = 


a>  —  1 


which  is  identical  to  Equation  (17)  of  Reference  3.  Also,  Equation 
(A-3)  becomes 


or 


03 


1  -  A 


a  +  ^  -  (i  -  a,  )  d  -  ttT)  } 


which  is  identical  to  Equation  (16)  of  Reference  3. 


Equations  (34)  and  (35)  of  Reference  1  may  be 


written : 


-2  .  ,,  x  x  rfXB 


(1  -  A^oT  -  2u>  +  (1  -  I2)  -  (o)  -  1)  =2 


=  0 


03  -  1  03  -1 

r  r 


(1  -  A1)  oj2  -  2oj  +  (1  -  I2)  +  (o32  -  1)  x 


XB(Uir  -  1)  _ 


=  0 


These  equations  can  be  rearranged  to  be  exactly  Equations  (A- 3) 
and  (A- 4)  above. 


■Thus,  it  is  seen  that  for  the  special  case  of  the 
isotropic  single-degree-of-freedom  hub,  the  equations  reduce 
to  the  standard  forms  obtained  by  others. 


